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ABSTRACT 

A thermocapillary parallel flow is established in a fluid filled 
slot with an applied temperature gradient. Low frequency jit- 
ter is imposed in arbitraiy directions. Vertical jitter proves 
to be relatively uninteresting, merely augmenting or opposing 
the basic thermocapillary flow. Streamwise jitter still pro- 
duces parallel flows, but these now exhibit boundary layers 
or layered cellular structures for large Rayleigh number as 
the applied stratification alternates between stable and unsta- 
ble. Runaways are possible for unstable stratification and 
these correspond to resonant excitation of stationary long wave 
Rayleigh-Benard modes. Spanwise jitter produces fully three 
dimensional motion. A spanwise-streamwise circulation re- 
sults for weak spanwise jitter, which advccts the interfacial 
temperature establishing a subsidiary spanwise thermocapil- 
lary flow. This flow is strong at small Biot number when 
advccted temperature is trapped in the slot, and has a counter- 
intuitive dependence on the spanwise-streamwise aspect ratio. 

1 INTRODUCTION 

Materials processing situations frequently involve interfaces 
with lengthwise temperature gradients, and are therefore sus- 
ceptible to thermocapillary motions. If materials processing 
is conducted in space, the thermocapillary flows will also be 
subject to fluctuating accelerations (p-jitter) characteristic of 
the space based microgravity environment [1, 2, 3, 4j. While 
there has been a large number of studies of buoyant systems 
subject to gravitational modulation [5, 6, 7, 8, 9, 10], there 
is no equivalent body of results for thermocapillary systems 
under these circumstances. 

The present paper considers the effect of jitter on a basic 
thcrmocapillary flow in a particularly simple geometry: a fluid 
layer of finite depth but infinite horizontal extent. This is the 
thcrmocapillary slot flow of Davis and coworkers [11, 12], 
which is sketched in Figure 1 . The underlying simplicity of 
the slot model is such that, even when jitter is added, the 
problem remains amenable to analytical techniques. This per- 
mits a thorough parametric study of the system, allowing the 
possibility of considering several different directions of jitter, 
including spanwise jitter which has been little studied to date 
even for buoyancy driven flows. The structure of the paper 
is as follows. In the next section we introduce the governing 
equations of the slot. In §3 we present results for jitter confined 
to the plane, and in §4 we consider spanwise jitter. Conclu- 
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Figure 1: Thermocapillary slot return flow produced by a 
lengthwise applied temperature gradient. A vertically varying 
advected temperature profile, which balances streamwise heat 
convection and vertical heat conduction, is also present. 

sions are given in §5. Full details may be found elsewhere 
[13,14]. 

2 GOVERNING EQUATIONS 

The governing continuity, momentum and thermal equations 
of the slot subject to g -jitter are 

V.u = 0, (1) 

MaPr~ l (u tt -F u.Vu) = - Vp + V 2 u -f- GT, (2) 

Ma(T, t + u.VT) = V 2 T, (3) 

where Ma is the Marangoni number, Pt is the Prandtl num- 
ber and G represents the dimensionless jitter, consisting of a 
dynamic Bond number G, multiplied by an 0(1) vector de- 
scribing the direction and time modulation of gravity. Note that 
the momentum equation can be recast as a vorticity equation 

Jie(u>,t -|- u.Vu;) = .Reuf.Vu + V 2 u> + VT x G, (4) 

which clearly shows the role of temperature gradients in vor- 
ticity generation. The temperature field is most conveniently 
written T = -x 4- T a (z,y,z,f) where -x represents the 
applied temperature field, that is responsible for the basic ther- 
mocapillary flow, and T a is produced by advection. 

The base of the slot y = 0 is an inpenetrable, no slip, 
thermally insulated boundary. The interface y — 1 is assumed 
to remain flat, to sustain thermocapillary stresses, and to permit 
heal transfer characterized by a Biot number Bi, Accordingly 

u = T y = 0 at y = 0 and 

v - u,y - 1 + T® = w <y + T* 

= T a y +BiT a = 0 at y = 1 . (5) 

The slot flows must also be consistent with no fluid penetration 
conditions which would apply on the sidewalls of a cavity 
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3 JITTER WITHIN THE PLANE 


of finite extent. There are various ways of ensuring this, 
depending on the class of slot flow. In many cases return flow 
conditions will apply 

j udy — j wdy~ 0. (6) 

J o J o 

However there are also more general possibilities (to be dis- 
cussed in §4). 

In the absence of jitter G = 0 the equations admit parallel 
flow solutions in the x direction, with the ad vec ted temperature 
field varying only in the vertical. The solutions are [II, 12] 

u=\y 2 -\ y, T a = -^Ma(3j, 4 -4j, 3 + l). (7) 

Our goal is to calculate how these solutions are modified in the 
presence of jitter. For simplicity we assume throughout that the 
modulation occurs at some frequency fi which is small enough 
to be quasistatic, allowing the neglect of time derivatives in 
equations (2)-(3). 

3 JITTER WITHIN THE PLANE 

If G is non-zero but confined to the xy plane, we still ob- 
tain parallel flow solutions, with advected temperature varying 
only in the vertical. It is clear from equation (4) that vertical 
components of gravity only produce vorticity when coupled 
to the applied temperature field, not when coupled to T a . As 
the applied temperature is independent of any thermocapillary 
motion, thermocapillarity and vertical gravity operate inde- 
pendently of one another. If G = (0, Geos fit, 0) vertical 
gravity produces fields 

u = -jjGcosQt (4 t/ 3 — ^y 1 4- 3y) , 

T“ = £MaGcosQ<( 5 V-iy 4 +iV 3 -|). (8) 

additional to those of equation (7). Thus vertical gravity alter- 
nately cooperates or competes with the thermocapillary fields 
as G cos fit executes its oscillation cycle. 

Now consider the case of streamwise gravity given by 
G = (Geos fit, 0, 0). From equation (4) we deduce that this 
generates vorticity when coupled with the advected tempera- 
ture field, implying an interaction will occur between thermo- 
capillarity and streamwise gravity. We define a streamfunction 
*0 such that u — . The thermal equation (3) now relates the 

streamwise heat advection Ma ij> and the vertical heat conduc- 
tion — T£. Substituting into vorticity equation (4), we deduce 

i>, vvvv — Ra cos fit = 0, (9) 

where the Rayleigh number Ra equals MaG. The boundary 
conditions are tp = = 0 at y = 0, and ip = ^ \ — 0 

at y = 1. 



Figure 2: Velocity profiles for streamwise jitter during the (a) 
stable and (b) unstable part of the oscillation cycle, with solid 
line: Ra\ cos fit| = 0, the base state profile; long dashes: 
Jla|cosfit| = 100; short dashes: J2a|cosfit] = 200; and 
dash-dot line: i?a| cosfit| = 2200. As Ra | cosfit| increases, 
the stable case (a) tends to form an interfacia] boundary layer, 
while the unstable case (b) forms cellular structures. Addi- 
tionally in (b), the profiles obtain a large amplitude if they lie 
near one of the runaway values of Ra cos fit. 

The solutions of this equation will be linear combinations 
of trigonometric and exponential functions. However these 
solutions are structurally quite different depending on the sign 
of Ra cos fit. If Ra cos fit is large and negative, the solutions 
tend to be confined in boundary layers (see Figure 2(a)). If 
Ra cos fit is large and positive, the solutions are highly cel- 
lular in nature (see Figure 2(b)). In both cases the thickness 
of the boundary layers or the wavelength of the cells scales 
as (Ra|cosfit|) 1/4 . There are also found to be particu- 
lar positive values of Ra cos fit which lead to infinitely large 
responses or runaways. 

Physical insights into the reason for the different behaviour 
in the two cases can be deduced by considering the strat- 
ified slot initially in the absence of thermocapillarity. If 
Ra cos fit < 0, gravity and the applied temperature gradi- 
ent are anti-parallel, meaning the slot is stably stratified in the 
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Figure 3: (a) Purely buoyant circulation flow induced by span- 
wise jitter in a slot with an applied temperature gradient. The 
circulation alternates between anticlockwise and clockwise as 
the jitter proceeds, (b) The circulation convects heat and ro- 
tates the isotherms (shown as dotted lines). The presence of 
a non-zero spanwise component of the interfacial temperature 
gradient drives a subsidiary thermocapillary flow. 


Rayleigh-Benard sense. If Ra cos fit > 0, gravity and the ap- 
plied temperature gradient become parallel, corresponding to 
unstable stratification. If thermocapillarity is now added to the 
stably or unstably stratified system, it breaks the symmetry of 
the Rayleigh-Benard basic state, thereby imposing directional- 
ity on the flow. The runaways in the unstable case correspond 
to values of the quasistatic Rayleigh number .Ra cos fit that 
support stationary, long wave Rayleigh-Benard modes, which 
become resonantly excited by the thermocapillary forcing. 


4 SPANWISE JITTER 

IfG = (0,0, G cos fit) equations ( 1 )— (5) no longer admit par- 
allel flow solutions, and in fact a non-linear set of equations 
must be solved and the flows are fully 3-dimcnsional. For weak 
spanwise jitter, we pursue a perturbation expansion in powers 
of Geos fit. It is clear from equation (4) that spanwise jit- 
ter coupling to the applied temperature field generates vertical 
vorticity at leading order. This produces a purely buoyant cir- 
culation flow around the slot (see Figure 3(a)). We concentrate 
solely on the leading order results here. Results through sec- 
ond order, which necessarily include steady streaming flows, 
are given in [14]. 


In determining this circulation, we want to continue to 
exploit the horizontal spatial invariance properties that follow 
from the geometric simplicity of the slot model. Accord- 
ingly we assume a first order circulation flow of the form 
G cos fit (ztioi(y),Oj 2^10(2/)) where «oi and idio are func- 
tions to be determined. These functions are necessarily non- 
return flows, and in place of equation (6) we do the following. 
Imagine that the infinite slot models a true shallow cavity of 
unit depth, length A~ l and breadth A ZX A~\ Consider a quad- 
rant 0 < x < |A“\0 < z < \ A zx A ~\ located between the 
cavity centre and one comer. No penetration conditions on the 
cavity sidewalls require that the fluid flux entering the quad- 
rant along matches that leaving along z — 0. Applied 

to the assumed slot flows this requires 


I w\ody = — A 2 zx I UQ\dy . (10) 

Jo Jo 

The parameter A has cancelled from equation (10) but the 
spanwise-streamwise aspect ratio A zx must be retained. The 
solutions for uq\ and u>oi are found to be 

“ 01 = -TTIirOv 2 -!'). ( ]1 > 

-.0 = T^r(W-v)- 02 ) 

Clearly streamwisc (spanwise) circulation flow is favoured at 
small (large) A tx . 

The circulation flow convects heat and therefore rotates 
the isotherms. This produces a temperature G cos fit zT 0 i (y) 
where 


Ma 

1 -F A 2 zi 


(hv 


W + 1 + 


skk 


(13) 


The term involving Bi accounts for the interfacial value of 
Toi and tends to be dominant even for Bi values as large 
as unity. A non-zero interfacial advected temperature arises 
here because the circulations are not return flows, and hence 
convect heat along the layer. If Bi — > 0 the convected heat has 
great difficulty escaping the interface, and so the temperature 
in the slot rises to a very high level. 

As there are now additional temperature gradients along 
the surface in excess of the applied temperature gradient, new 
subsidiary thermocapillary motions are produced (see Fig- 
ure 3(b)). These will be in the spanwise direction and we 
write them as G cos fit u>oo(t/), where 


Ma 


U700 — 


(1 +A\ I )Bi 


• (s» - b 1 ) ■ 


(14) 


If thermocapillarity is inherently strong (Ma is large) or if 
heat is trapped in the slot (Bi is small) the subsidiary flow 
can dominate the circulation. The parametric dependence of 
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woo on A zx is counterintuitive on geometric grounds, i.e. the 
span wise subsidiary flow is largest when spanwise-streamwise 
A zx aspect ratio is small. The reason is that small A zx favours 
the streamwise circulation that produces the isotherm rotation 
which in turn drives woo. 

5 CONCLUSIONS 

The thermocapillary slot model subject to ^-jitter exhibits rich 
physical behaviour. The simplicity of the slot geometry allows 
an analytic approach to the problem, facilitating a thorough 
parametric investigation, including the possibility of consider- 
ing various directions of jitter. Vertical jitter is uninteresting 
because the advected temperature fields that arc produced do 
not themselves produce any vorticity. Streamwise jitter is more 
interesting, and produces velocity profiles exhibiting bound- 
ary layer structures in the stable part of the cycle and cells in 
the unstable part. In the latter case it is possible to resonantly 
excite stationary long wave Rayleigh-Bdnard modes for par- 
ticular values of the quasistatic Rayleigh number. Spanwise 
jitter generates circulation flows around the slot. As the circu- 
lation advccts interfacial temperature, it establishes subsidiary 
spanwise thcrmocapillary flows. At large Marangoni number 
or small Biol number, these subsidiary motions can become 
important. 
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